Copolymer Networks: 
Multifractal dimension spectra in polymer field theory 

C. von Ferber 1 ' 2 and Yu. Holovatch 3 
1 Fachbereich Physik, Universitdt - Gesamthochschule - Essen, 
D-45117 Essen, Germany 
2 School of Physics and Astronomy, Tel Aviv University, IL-69978 Tel-Aviv, Israel 
3 Institute for Condensed Matter Physics, Ukrainian Academy of Sciences, 
t> ' U A- 2900 11 Lviv, Ukraine 

os : 

, We explore the rich scaling behavior of copolymer networks in solution. 

We establish a field theoretic description in terms of composite operators. Our 
3rd order resummation of the spectrum of scaling dimensions brings about 
remarkable features: Convexity of the spectra allows for a multifractal inter- 
pretation. This has not been conceived for power of field operators of (f) A field 



(N 



theory before. The 2D limit of the mutually avoiding walk star apparently 
corresponds to results of a conformal Kac series. Such a classification seems 
not possible for the 2D limit of other copolymer stars. The 3rd order calcu- 
^C} ■ lation of a large collection of exponents furthermore allows for a consistency 

check of two complementary schemes: epsilon expansion and renormalization 
\ at fixed dimension. 

o 
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I. INTRODUCTION 



e 

Recently much interest focused on the relation of field theory and multifractals HB 
and the associated multifractal dimension spectra [^,|| as well as non-intersecting random 
walks and their 2D conformal theory ||. We present a model of multicomponent polymer 
networks that shows a common core of these topics and allows for a detailed study of the 
interrelations. The flux of diffusion onto an absorbing fractal defines a multifractal measure. 
Cates and Witten have mapped the moments of this flux to that of a star of random 
walks (RW) avoiding the absorber taken to be a polymer or RW itself. Using the field 
theoretic formulation of polymer theory we show that the spectrum of scaling exponents 
governing these problems is given by the anomalous dimensions of composite operators with 
appropriate symmetry. 

For polymer networks consisting of polymer chains of one species it has been shown, that 
the basic scaling exponents are connected with 'stars', polymer chains tied together at one 
core The number of configurations Z*f of a polymer star with / arms of TV monomers 

will scale for large N like 

z* f ~ ah''- 1 ~ {Riiyi-fo. (1) 

The second part shows scaling with the size R ~ N u of the isolated coil of N monomers on 
some scale I. The exponents v = 3/4,0.58(8) and 71 = 72 = 7 = 43/32, 1.16(0) for space 
dimensions d — 2, 3 are known in polymer theory |J. The exponents 7/ have been calculated 
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analytically in perturbation theory []7|,|8|,p!0[| , by exact methods in two dimensions |J, and by 



Monte Carlo simulations 111 



At short distance two polymer stars will repel each other. In view of the below advocated 
language of field theory this is described in terms of a short distance expansion. One finds 
the following relation for the probability P(r) to find the cores of two stars of f\ and f^ at 
short distance r @ 

P{r) ~ r e , 6 = Vfl + Vh - Vfl+h > . (2) 

This is compatible with the result, that the spectrum of polymer star exponents rjf is convex 
from below as function of / with rji = 0. 

On the other hand a multifractal (MF) measure /i x defined on the sites x of scale I on 
some object of size R is characterized by the scaling of its moments averaged over all sites: 

<a£> = £a£~(W- (3) 

X 

From general inequalities for the moments of a probabitity distribution one may deduce 
that the spectrum of exponents yf has to be convex from above. This indicates an apparent 
discrepancy between objects described in field theory (FT) as powers of field (see below) 
such as polymer stars, and the moments of a MF measure 0. This we want to resolve by 
including both concepts in the same FT formalism showing that they are special cases of 
a more general approach, which in addition also describes the problem of non-intersecting 
random walks. 

To this end we study the scaling behavior of a polymer star or a general network of chains 
of different species and thus, within a unique formalism, include effects caused by self and 
mutual interactions between polymers of different species forming a network. We combine 
the field theoretic formalism developed for the description of polymer stars and networks || 
with the corresponding theory which describes multicomponent polymer solutions W\ . 



II. THEORY 

We introduce a Landau-Ginsburg-Wilson-Lagrangian C of / interacting fields <f>\, each 
with n components, i.e. (f) 2 a = Sa=i(0a) 2 5 with an interaction matrix u aa i and mass param- 
eters m a \ 

f f 
C{<P b ,m b } = - £ jd d r (m a <f>l + (V0 a (r)) 2 ) + - £ u aa , J d'r#(r)#(r). (4) 

a=1 ' a,a'=l 

In this theory the star exponents are given in terms of the anomalous dimensions of composite 
operators ]l{=i <fia ||- We define vertex functions T*f with insertion of this operator by 

5(q + ... + q f )T*f(q . . . q f ) = f {[ e^d d r k ({[ U^Un) . . . ^/(r,))^, (5) 

fc=0 a=l 

As in standard polymer FT this is evaluated with respect to the Lagrangian (f|) keeping 
only contributions which correspond to one particle irreducible (lpi) graphs which have 
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nonvanishing tensor factors in the n = limit. In the single component case the theory may 
also be described in terms of one 0(n) symmetric field with n > f, where the corresponding 
operator is N ai '" a f (f) ai ■ ■ ■ <p af with a traceless tensor N ai --- a f in the formal limit n = P,|l3"f. 

We apply RG theory to make use of the scaling symmetry of the systems in the asymp- 
totic limit to extract the universal content and at the same time remove divergences which 



occur for the evaluation of the bare functions in this limit ||14j| . Several asymptotically equiv- 
alent procedures serve to the purpose of renormalization. In the present study we use two 
somewhat complementary approaches: zero mass renormalization with successive e = 4 — d 
-expansion |TJ]] and the massive RG approach at fixed dimension [T5| . Application of both 
approaches will enable us to check the consistency of approximations and the accuracy of 
the results obtained. We pass from the theory in terms of the initial bare variables to a 
renormalized theory. This can be achieved by a controlled rearrangement of the series for 
the vertex functions @ introducing renormalizing Z-factors for fields (Z^), couplings [Z ab ] 
and mass. Then, for instance the bare couplings u ab are given in terms of their renormalized 
dimensionless counterparts g ab by 

u a b = ft 4 d Z(j> a Z ( f >b Z ab g ab . (6) 

The scale parameter k represents the mass at which the massive scheme is evaluated and the 
scale of external momenta in the massless e-expansion scheme. We define the Z-factors in 
(§) as to renormalize the correlators (■ • •) in each RG procedure (see e.g. ]nj). The polymer 
limit n = of zero component fields leads to essential simplification. Each field <p a , mass 
m a and coupling u aa renormalizes as if the other fields were absent. The renormalization of 
the couplings u ab involves only the fields </> a ,</>6 fl^| . The renormalized couplings g ab defined 
by relations (0) depend on the scale parameter k. Thus the renormalization Z - factors 
also depend implicitly on k. This dependence defines the RG functions and exponents: 
K-^9aa = Paa(gaa); Kgjjflto = Pab(9aa, 9bb, 9ab); K-^lnZ^ = The function rfa a defines 

the pair correlation critical exponent. The set of scaling exponents 77*/ for general copolymer 
stars is defined by the renormalization factors Z*f for the star vertex functions T*?: 

f d 
Jl Zj^Z^ f Y* f {u w {g bb ,g Vbl ,g w )) = k s s , with r] iff (g ab ) = K—lnZ^ . (7) 
0=1 aK 

5f — d+ (1 — d/2)f is the engineering dimension of the corresponding bare vertex function. 

In a study devoted to ternary polymer solutions the RG flow given by the above defined 
/3-functions has been calculated flr2| , |TE| to third loop order. The equations for the fixed 
points of the /^-functions were found to have the following nontrivial solutions: /3 aa {9s) = 

and for a ^ b: ^(0,0, g* G ) = 0, 0, 0u) = 0, MO, 9* s , 9*v) = 0, 0s) = 0, 

corresponding to all combinations of interacting and non-interacting chains. 

We evaluate the exponents for two general arrangements of the fixed point matrix. The 
ternary case of two mutually interacting species of polymer chains in solution, and the 
mutual avoiding walk case of essentially / only mutually interacting species. In the first case 
we describe polymer stars made of fi chains of species 1 and / 2 = / — fi chains of species 
2. Either both species are non self-interacting and 

^/i/ 2 = V*f(9ab = if a, b < }\ or a, b > fi, else g ab = g£) , (8) 
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or species 1 self-interacts and species 2 does not such that 

"^hh = VtfiQab = 9s if a, b < /i; g ab = if a, b > f X ] else g ab = gfa). (9) 

For $i = this includes the homo-polymer star with rjf — rfj Q in eq.(|l]). The mutually 
avoiding walk case reads 

r]f AW = = if a = b else g ab = g* G ) . (10) 



III. RESULTS 

We give the results for the exponents in e = 4 — ci-expansion. The corresponding more 
lengthy expressions obtained by fixed d = 3 RG may be found in |Tj| : 

Vhf 2 (z) = Si h\+h U {h - 3 +/,) j /* (f s - 3 + /,)(/, +h + 3C(3) -3)^ (11) 



n u 



u hh {e) = /, (l - U - 3/ 2 ) ~ + /, (25 - 33/, + 8/, 2 - 91/* + 42/; f 2 + 18 f 2 2 )^ 

-/, (577 - 969 /, + 456 /, 2 - 64 /, 3 - 2463 f 2 + 2290 /, / 2 - 492 /, 2 f 2 + 1050 f 2 2 
-504/, // - 108// - 712 C(3) + 936/, C(3) - 224/, 2 C(3) 

+2652 f 2 ((3) - 1188/, /,C(3) - 540 /, 2 C(3))^ (12) 

r]f m (e) = -(/ - l)/ \ +f (/ - l)(2/ - 5)^ - (/ - l)/ (4/ 2 - 20/ + 8/C(3) - 19C(3) + 25)^ 

(13) 

Here ((3) ~ 1.202 is the Riemann ^-function. The above formulas reproduce the 3rd order 
calculations of 7/ — 1 = ^(i]^ — fr^o) f| as well as the 2nd order exponents A'**) defined in 
equations (xx) of §, \^(n) = - % G n , \^(n) = -% u n + % u , A( 48 )(n) = -7^, X^(n) = 
—rjYn, correcting a missprint in eq.(49) of |||. Also the 2nd order results for exponents 
xl,u — = —^(VLn ~ Vli) °f [0 an d a L = 1/277^ AW defined in |J find their 3rd order 
extension by the above expansions. 

With these exponents we can describe the scaling behavior of polymer stars and networks 
of two components, generalizing the relation for single component networks @. In the 
notation of ([I]) we find for the number of configurations of a network Q of F\ and Fi chains 
of species 1 and 2 

Z g ^(R/e)^- F ^- F ^\withr ]g = -dL+ £ NhhVhk, (14) 

/1+/2M 

where L is the number of Loops and N^f 2 the number of vertices with f\ and f'2 arms of 
species 1 and 2 in the network Q. To receive an appropriate scaling law we assume the 
network to be built of chains which for both species will have a coil radius R when isolated. 
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To obtain reliable numerical values from the e-expansions in ([11]) - (|13[) and from the 
series obtained in the fixed d scheme [[l(| we apply Borel resummation using the technique of 
conformal mapping [T7|] which has proven to yield good results for many critical exponents. 
We use information about the higher order behavior of the series (|TT|)-(|13|) derived 

from the instanton analysis of the appropriate field theory. The results for d = 3 are given 
in Table |. The data show consistency and stability of the results while deviations grow for 
large number of arms as may be expected. Note that the above expansions are in fact series 
in fe, not e alone. 



IV. CONCLUSIONS 

A. Multifractals and Field Theory 

Does the data answer the question of convexity? A close study of the matrix of values 
reveals, that for fixed fi both rff^ 2 and are convex from above as function of f%, thus 
yielding 'MF statistics'. The relation to a MF spectral function for j\ = 1,2 has been 
pointed out in |J, it is analysed in close detail in view of the new data and FT formulation 
in a separate publication fl6| . On the other hand also copolymer stars should repel each 
other. This is found to be true as well, the corresponding convexity from below shows up e.g. 
along the diagonal values r/ff as function of /. The general relation r nh}2 Jrr lf 1 f 2 — Vfi+f'Ju+fz 
is always fulfilled. In view of our FT formalism the MF moments (fi k ) are represented by 
field operators (fi^b = 0<u ' ' ' 0a L 0i>i • • • <f>b k i n a FT with vanishing interactions g^b ■ Thus, 
even though simple power k of field operators 4> k do not describe MF moments (T| , they may 
be written as a power L + k of field operators which have the appropriate short distance 
behavior. This is also illustrated in fig.l, showing the spectrum of exponents in the 2D 



limit [16]. The opposite convexity along the two axes is clearly seen for these unsymmetric 



combinations of a polymer /x-star and a random walk / 2 -star which mutually interact. 



B. 2D Copolymer Stars 

The 2D exponents for polymer stars have been shown to belong to a Kac series of 
exponents of conformal FT with 7^ — 1 = (4 + 27 / — 9/ 2 )/64 ||. There are strong 
indications that this is the case also for MAW stars with ^¥ AW = (1 - 4/ 2 )/12 0. 
Already in view of fig.l though, such a simple 2nd order polynomial seems not to de- 
scribe the 2D limit of general copolymer star exponents. In 2D however, each chain 
of a star will interact only with its direct neighbors. A star described here by rfj^ 
will behave like a MAW 2/-star if each species-1 chain has two neighbors of species- 
2 whereas it will behave differently if the chains are ordered such that each species is 
in one bulk of chains. The 2D copolymer stars in this sense reveal an even richer be- 
havior. Thus, the copolymer generalization of the MAW star adds another problem, for 

which a rigorous formulation in terms of an exactly solvable 2D model is yet to be found. 

*** 

We thank Lothar Schafer for valuable discussions. Supported in part by Deutsche 
Forschungsgemeinschaft SFB 237 and Minerva Gesellschaft. 
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TABLE I. Values of the copolymer star exponent ffj x ^ upper part (U), and ??^j 2 > lower part 
(G), at d = 3 obtained by e-expansion (e) and by fixed dimension technique (3d). 
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1 


-0.43 


-0.45 


-0.79 


-0.81 


-1.09 


-1.09 


-1.35 


-1.37 


-1.60 


-1.64 


-1.81 


-1.89 




2 


-0.98 


-0.98 


-1.58 


-1.60 


-2.13 


-2.19 


-2.61 


-2.71 


-3.05 


-3.21 


-3.46 


-3.68 


u 


3 


-1.64 


-1.67 


-2.44 


-2.52 


-3.16 


-3.30 


-3.82 


-4.04 


-4.44 


-4.75 


-5.01 


-5.42 




4 


-2.39 


-2.47 


-3.33 


-3.50 


-4.20 


-4.48 


-5.02 


-5.40 


-5.80 


-6.30 


-6.53 


-7.15 
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-3.21 


-3.38 


-4.28 


-4.57 


-5.28 


-5.71 


-6.24 


-6.81 


-7.15 


-7.89 


-8.02 


-8.92 
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-1.35 


-1.63 


-1.69 


-1.88 


-1.98 


-2.10 


-2.24 
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-2.53 
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-3.17 


-3.51 


-3.75 


-3.95 


-4.28 
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-4.50 


-4.94 


-5.36 


-5.62 


-6.15 
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-5.27 


-5.71 


-6.24 


-6.84 


-7.12 


-7.90 
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-7.42 


-8.24 


-8.50 


-9.54 
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-9.78 


-11.07 




FIG. 1. Exponent f?}^ in the 'Unsymmetric' fixed point at d = 2 obtained in e-expansion and 
in fixed d scheme. The steps in the 'flying carpet' indicate the difference of the results in the two 
approaches 



7 



